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TECHNICAL PAPER 


A STUDY TO EVALUATE STS HEADS-UP ASCENT TRAJECTORY 
PERFORMANCE EMPLOYING A MINIMUM-HAMILTONIAN 
OPTIMIZATION STRATEGY 


I. INTRODUCTION 


The Space Transportation System (STS), more commonly known as the Space Shuttle, is currently 
America’s primary heavy-lift launch vehicle. Therefore, any vehicle performance enhancements are of 
significant interest. One method identified to potentially increase STS performance is to alter the trajectory 
of the vehicle such that it flys in a “heads-up” mode. The term heads-up refers to the relative position of a 
crew member in the shuttle with respect to the ground (Fig. 1). 

The objective of this study is to ascertain the performance capability of the STS, as a function of 
Solid Rocket Motor (SRM) bumrate and maximum dynamic pressure (max Q), for heads-up flight. The 
results can then be compared to the current heads-down flight mode. 


II. METHOD 


A. Simulation Program 

The program used to simulate shuttle ascent trajectories, in this study, is referred to as the Minimum 
Hamiltonian Ascent Shuttle Trajectory Evaluation (MASTRE) program. Development of this program was 
completed in November 1974 as an outgrowth of three NASA Marshall Space Flight Center (MSFC) 
programs: RAGMOP [1], ROBOT [2], and Lifting ROBOT [3]. Trajectory optimization is obtained by a 
steepest ascent algorithm using a minimum Hamiltonian (min-H) strategy, where the derivation has been 
expanded to encompass the effects of the atmosphere and of moment balance control. The combined effects 
of aerodynamics and fixed engine offsets are counterbalanced, or moment balanced, by controllable 
engines. Three degrees-of-freedom are maintained in the program by instantaneous command and response 
of the controllable engines; therefore, no vehicular angular rates are considered. Note, the controllable 
engines can be either the SRMs or the Space Shuttle Main Engines (SSMEs), but not both. 

The pitch attitude angle (\p), yaw attitude angle (xy)> and roll attitude angle (xr) used in the 
MASTRE program are defined in Figure 2. The X P , Y P , and Z P axes form an Earth-centered inertial 
plumbline coordinate system; where Y P is the local vertical, X P is directed along the flight azimuth, and Z P 
forms a right-handed system. The body coordinate system is defined by X m , Y m , and Z m axes; where Y m is 
along the longitudinal axis of the vehicle, X m points away from the orbiter tail, and Z m forms a right-hand 
set. Therefore, x P is measured from the local vertical to the projection of the external tank centerline in the 
Y P -X P plane. Whereas, xy is measured from the ET centerline to the projection of the ET centerline in the 
Y P -X P plane. Finally, xr is measured about the ET centerline. 




) 


y p 



Figure 2. Attitude angle definition. 

The gravitational model in the program is derived from the rotating Fischer Earth model, which 
establishes a “pear shaped,” ellipsoidal Earth model. The model also serves as the reference for altitude 
computation. Details of this model are described in Reference 4. 

The other models and/or groundrules used in this study are listed in the Groundrules section of this 
report. A description of the min-H optimization process used by the MASTRE program is presented in the 
next section. 
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B. Optimization Methodology 


The trajectory optimization procedure employed in the MASTRE program is known as the 
Minimum Hamiltonian Method. A basic discussion of this optimization scheme follows. 

Given, a set of n nonlinear differential equations 

X; = fj (u £ , X K , t) ; i, k = 1 — n ; £ = 1 — r (2.1) 


where 


u £ - control variable (not differentiated) 

X K - differentiated variable (one derivative only) 
t - time . 


In the MASTRE program, the nonlinear differential equations are the primary equations of motion. These 
equations are integrated using a fourth-order Runge-Kutta numerical integration scheme and are given by 
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(2.2) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


4 



Where, 


F - the total force on the vehicle due to thrust aerodynamics and winds 
M - the total vehicle mass 
G - the gravitational acceleration 


Note, equation (2.8) is not an equation of motion and does not require numerical methods for integration. 
The flow rates are input data and are shown to point out the source of vehicle mass changes. 

Also, 


SRM Thrusts - are obtained from user input sequential access files 

SSME Thrusts -are tag values, which are input via the ENGDAT namelist variable (see Appendix) 

Aerodynamic Forces - are determined by table look-up, as functions of Mach number, angle of 

attack (a), and angle of sideslip (P) 

Wind Forces - are determined by table look-up as a function of altitude, and the wind month is input 
by the namelist variable IWNUM (see Appendix). 

Now, rearranging the equations of motion, the following equation is established 


yi = Xj - f ( 


i = 1 


(2.9) 


Then, the introduction of Lagrange multipliers leads to 


E = \j(t)yi = \i(t) [Xj-fi] * i = 1 — * n 


(2.10) 


From calculus of variations, an integral can be minimized (or maximized) via the use of the Euler equation , 
which results in the following : 
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which leads to, 



S 


K= 1 




dXi 


= 0 


or, 


n 

v . y v 5f K . 

“h ,v — 0 * i 1 *n 

K = 1 ^ x i 


1 — ► n 


( 2 . 12 ) 


(2.13) 


and, in matrix form 


X = 



(2.14) 


This equation is known as the adjoint equation. 
Also, 


d dE dE 

dt du du 


(2.15) 


which yields, 


2 


K = 1 


du 0 


= 0 


where, 


/ = 1— r 


(2.16) 



(2.17) 


So, in matrix form 



(2.18) 
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or. 



(2.19) 


where, H = f'X is known as the Hamiltonian. The above equation is an /x 1 matrix of the partial deriva- 
tives of the Hamiltonian with respect to control variables, and is the optimality equation. MASTRE 
program control variables used in this study include Xp, Xy, main engine cutoff (MECO) time, and SSME 
throttle down to 65 percent initiation time during first stage. 

Initial conditions must be satisfied in the solution process, and these conditions are given by 


Xj(to) = X 0 , , i = 1 — n 


( 2 . 20 ) 


These variables include initial vehicle position, velocity, acceleration, and mass. 

Also, the P terminal or intermediate constraints have to be satisfied. These constraints are given by 


M [X(t),t] , j = l — p 


( 2 . 21 ) 


where. 


tytEco. for terminal constraints 
tj, for intermediate constraints 


The terminal constraints are MECO targeting conditions, which include altitude, velocity, and 
flight path angle. In the case of intermediate constraints, time is no longer MECO time but some 
intermediate time. Intermediate constraints used in this study are vehicle dynamic pressure and vehicle 
longitudinal acceleration. 

Now, while satisfying all of the previous conditions, a payoff quantity is to be minimized or 
maximized. This quantity is given by 


<J> — <}>[X(tMECo)*tMECo] 


( 2 . 22 ) 


In this study the payoff quantity is the MECO or injection weight, and this parameter requires 
maximization to achieve maximum performance. 
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Note, the Lagrange multipliers used in the adjoint and optimality equations are given by 


MO = MO + M^ v 


(2.23) 


where 

X - is a nxl matrix 

X , - is a nx 1 matrix 
9 

X^ - is a nxP matrix 

v - is a Px 1 matrix and consists of constant Lagrange multipliers associated with the terminal or 
intermediate constraints 



which is evaluated at t = t meco 



which is evaluated at 


(2.24) 


(2.25) 


1 tMECO> for terminal constraints 
t = J t[, for intermediate constraints 


The above equations are obtained by using the transversality conditions, and again they allow for the deter- 
mination of the Lagrange multipliers at the terminal and/or intermediate time. 


C. Optimization Iteration Process 

For optimum conditions, as shown in equation (2.19), H u = 0. However, in general this relation 
will not be true for initially chosen values of the control parameter, u(t), and the Lagrange multipliers 
associated with the terminal or intermediate constraints, v. Therefore, in order to satisfy the optimaiity 
equation (2. 1 9), initial conditions (2.20) and terminal/intermediate constraints (2.21 ), an iterative solution 
process is necessary. 
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The first step in this procedure involves the selection of initial values for the control parameters (\p, 
Xy, MECO time, and SSME throttle down to 65 percent initiation time during first stage) and the Lagrange 
multipliers associated with the terminal or intermediate constraints. The optimization method generates 
changes to these initial values, and this process is described by 


(i+D 


u(t) = 'u(t) + 8u(t) 


(2.26) 


and 


(i+D 

v 


'v + dv 


(2.27) 


where, 


i - is a superscript representing the i th iteration 

(i + 1 ) - is a superscript representing the (i lh + 1 ) iteration 

8u(t) and dv - are zero for the 0 lh iteration. 


Now, knowing the above values, a forward integration of the n nonlinear differential equations (2.1) 
occurs from X(t Q ) = X 0 to the currently established MECO time. Next, the adjoint equation (2.14) is 
integrated backwards from the known MECO time to t c , in order to calculate the Lagrange multipliers. 
Note, since the MECO time is known, the initial value of k is obtained from equation (2.23), where t = 
tMEco Also, during each step of the backward integration of the adjoint equation, the value of the opti- 
mality equation is calculated from 


H u 



(2.28) 


which is a combination of equations (2.18) and (2.19), and for optimum conditions H u = 0. Recall, in 
general the above equation will not be zero for the initially chosen u(t) and v. Therefore, to force H u to zero, 
changes need to be made to the control parameters, u(t), and the effect these changes have on the constraints 
need to be accounted for by also changing the Lagrange multipliers associated with the constraints, v. 
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/ *M 

fa 


X j, G(u* - u c ) dt 


and 


l dP 2 - dv];' 14 dvA 1/2 

l Icc — ^>l) c I^4c J 


Now, 


J l MECO 

, 5u' H uu 5u dt 

to 


and 


C l MECO 

Icc = Ju, (u* - Uc)' H uu (u* - u c ) dt 


(2.33) 


(2.34) 


(2.35) 


(2.36) 


The strategy involved in determining 5u(t) and dv is to first select a value for dP 2 and set dvp = -i|>. Now, for 


dP 2 -df i* diV < 0 
then set P = 0 and dvj> = -Ki|>, where 
dP 2 

d*’ I** di|> 

or for 

dP 2 - dvV' I^dili^O 

then set P = 1 or to the calculated P, whichever is the smaller value. Once P = 1 , discontinue the use of dP 2 
to calculate P. 

Hence, with 5u(t) and dv known, the iteration procedure is repeated until the optimality equation, 
initial conditions, and the constraints are satisfied within some tolerance. 




The process just described is again referred to as chord min-H and is a convergence scheme 
designed for use with high-speed digital computers. Further information on the chord min-H formulation is 
described in Reference 5. 


D. MASTRE Program Modifications 

As stated in the introduction, the purpose of this study is to examine the performance implications 
of heads-up shuttle ascent. However, previous trade studies conducted by Rockwell International [6] 
indicate that rolling the vehicle to heads-down, from an initially heads-up ascent, prior to Solid Rocket 
Booster (SRB) separation increases STS performance. This performance increase occurs because if a 
heads-up profile is maintained during staging, thrust from the SSMEs (which provide trim control during 
SRM thrust tailoff and after staging) is forced to a non-optimum angle. The non-optimum condition is a 
result of angle of attack (a) constraints on the vehicle and cause roughly a 900 lb performance loss. Also, a 
heads-down attitude after SRB separation is preferred for intact abort considerations. So, two vehicle roll 
maneuvers are required for what is referred to as a heads-up trajectory in this study. The first roll maneuver 
occurs after the launch pad is cleared to get the vehicle in a heads-up position; then a second roll maneuver 
is performed prior to SRB jettison to place the vehicle in a heads-down attitude. Originally, only the single 
roll maneuver after launch pad clearance was allowed to occur in the MASTRE program. Therefore, an 
additional roll option needed to be incorporated into the program, and as a result two new inputs were 
created as follows: 

1 . HDMACH - Mach number to initiate second roll maneuver. 

2. CHRDHD - roll rate for the second roll maneuver in degrees per second. 

Note, the second roll maneuver is set up in the program to go from a heads-up attitude to a heads-down 
attitude only. Otherwise stated, a heads-up attitude to a heads-sideways attitude cannot occur; the vehicle 
will always roll to heads-down. 

Another modification to the program involves the calculation of SSME specific impulse (Isp). The 
Isp’s of the SSMEs is not considered constant but is assumed to vary as a quadratic function of engine rated 
power level (RPL) as follows: 


SSME Isp = C, P L 2 + C 2 P L + C 3 


(2.38) 


where 


C|,C 2 ,C 3 - are coefficients of constant value 
Pl - engine power level 
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The previous version of the program required the user to input the coefficients of the above equa- 
tion. However, the program has been modified to calculate the coefficients based on three input Isp’s and 
their respective power levels. The input variable names are TBLISP and TBLPL for the SSME Isp’s and 
power levels, respectively. This procedure is useful because SSME assessment tag Isp’s are given at three 
different power levels. Overall, this program modification was used to simplify user input. 

A final minor modification to the MASTRE program dealt with the input of SRM ballistic perform- 
ance data (pressure, thrust, flowrate and nozzle throat exit area time histories). The old version of the 
program required the SRM performance data to be entered in a namelist format (see the Appendix as an 
example of namelist input). This input method has been changed so the MASTRE program can now read a 
sequential access file containing the SRM data, and only the name of the file requires input. Therefore, the 
process of running trajectories using a variety of SRM performance characteristics can easily be accom- 
plished, as required in this study. 


E. Trajectory Simulation Groundrules 

The groundrules used in an ascent trajectory simulation must be explicitly stated because of the 
multitude of parameters which can potentially affect the results. However, to include all the groundrules 
requires a great deal of detail. Therefore, the following is a list of documents and the groundrules they 
contain, which were necessary for this study. 

1 . Shuttle Operational Data Submittal (SODS) J-763 Addendum 1 [7] — This document provides 
all the data required for simulating performance reference mission four (PRM-4). This mission is a baseline 
trajectory, launched in a polar orbit from the Western Test Range (WTR), and is used for trajectory design. 
However, not all the inputs for this study were taken from this SODS. The reason for this decision is 
because some of the data is outdated, and SRM data had to include the effects of various propellant 
bumrates (Rt>’s) as described in item 5. The following information was taken from this SODS: 

a. Launch Site - WTR. 

b. Orbit Destination - 150 nautical miles (n.mi) altitude, inclined 98 deg to the equator. 

c. MECO Targets - Altitude of 57 n .mi . , relative velocity (V R ) of 25,374 ft per sec , and relative 
flight path angle (y R ) of 0.65 deg. 

d. STS Vehicle Aerodynamics - Integrated Vehicle Baseline Configuration 3 (IVBC-3) 
aerodynamics. 

e. Launch Site Winds - December mean monthly wind model for WTR. 

f. Launch Site Atmosphere - December Range Reference Atmosphere (RRA) for WTR. 

g. Vehicle Weight Summary Data (except SRM weights). 

h. Propellant Loads - Orbital Maneuvering System (OMS) and Reaction Control System 
(RCS) loads only. 
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i. Longitudinal Acceleration - Constrained to a maximum value of 3 G’s. 

2. Rockwell International Internal Letter, FSD&P/AFS-85-36 1 [8] — This letter contains inform- 
ation for a heads-down, August launch trajectory, using filament wound case (FWC) SRMs with a bumrate 
of 0.37 1 in. per second (ips). Note, all bumrates in this study are referenced to standard conditions of 60°F 
and 625 psia. Now, accompanying this letter was information for a heads-up, December launch trajectory, 
using FWC SRMs with 0.392 ips R b ’s. Data for this study was taken from this accompanying heads-up 
trajectory. Specifically, these data include angle of attack (a) and angle of sideslip ((3) Mach number his- 
tories, which are consistent with the nominal design requirements established in Reference 7. Also, the 
vehicle center of gravity (eg) Mach number history was taken from the heads-up trajectory. Note, this eg 
data incorporates FWC SRMs with R b ’s of 0.392 ips. 

3. Main Propulsion System (MPS) Propellant Inventory STS 61-1 Revision W [9] — This docu- 
ment contains the propellant loads for the ET. The propellants are liquid oxygen (LOX) and liquid 
hydrogen (LH 2 ), and are used by the SSMEs. 

4. MSFC Memo EL24 (86-36) [10] — Data contained in this memo are the MPS assessment tags. 
These data include SSME thrust, flowrate and Isp at various power levels. Note, SSME vacuum thrust at 
100 percent RPL is 469,760 lbf. 

5 . Heads-Up Vandenberg Launch Site ( VLS) FWC High Bum Rate SRM Data Submittal [11] and 
STS-2V(62-B) SRM Data Submittal [12] — These data submittals contain propellant load, mass property, 
and ballistic performance scaling information for FWC SRMs with R b ’s of 0.368 ips and 0.392 ips, respec- 
tively. The information in these data submittals was used to linearly interpolate/extrapolate all the SRM 
data required in this study. Note that to increase the bumrate of an SRM, iron oxide is added to the propel- 
lant. This procedure increases the propellant weight but reduces propellant Isp. Also, because an increased 
bumrate leads to increased motor chamber pressure, stronger and heavier motor cases are required. There- 
fore, SRM inert weight also increases with increasing bumrate. The specific linear equations used to 
account for the above three effects are as follows: 

a. Propellant Weight 

WP = 89208 R b + 1074388 (2 39) 

where, 

WP - is propellant weight in pound mass (lbm) 

R b - is propellant bumrate in ips. 

b. Inert Weight 

WI = 84542 R b + 129546 (2.40) 

where, 

WI - is SRM inert weight in lbm. 


14 



c. Isp Reduction Multiplier 

M = - 0.0810583 R b + 1.0298295 (2-41) 

where 

M - is the ISP reduction multiplier. 

Note, that even though the SRM mass properties change with bumrate, the STS vehicle eg Mach number 
history was assumed constant on all trajectory runs. This assumption was made because the overall vehicle 
eg is not significantly changed, and so the effect on STS performance is negligible. Again, the eg data for 
the vehicle was taken from the heads-up trajectory accompanying Reference 8. 

The SRM ballistic performance data (time histories of chamber pressure, thrust, flowrate, and 
nozzle throat exit area) were scaled from an SRM Block II baseline motor designated as TC-MD-27 1 -84- 
MAX2. This baseline motor had to be scaled to account for bumrate, propellant weight, and propellant 
mean bulk temperature (PMBT). These empherical scaling equations are as follows: 


THRUST = (THB) (e) Kl (KW) 1 53846 

(2.42) 

PRESSURE = (PB) (e) Kl (KW) 1 53846 

(2.43) 

TIME = (TIB) (e) K2 (KW)' 0 53846 

(2.44) 

FLOWRATE = (FLB) (e) K2 (KW) 1 53846 

(2.45) 


where THB, PB, TIB, and FLB are thrust, pressure, time, and flowrate of the baseline motor, respectively. 
Also: 


K, = 0.0011 (PMBTS - PMBTB) + 1.567153 In (RB/RBB) 

K 2 = 0.001063 [(PMBTS - PMBTB) + 1474.274 In (RB/RBB)] 

KW = WP/WPB 

where 

PMBTS - PMBT at the launch site in degrees Fahrenheit (depends on launch date) 


(2.46) 

(2.47) 

(2.48) 
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PMBTB - PMBT of the baseline motor in degrees Fahrenheit 

RB - motor bumrate desired in ips 

RBB - baseline motor bumrate in ips 

WP - propellant weight of motor desired in lbm 

WPB - propellant weight of the baseline motor in lbm. 

Once a motor was properly scaled, then the Isp reduction multiplier was applied. Also, SRB staging time 
was set to occur 6.2 sec after the time when chamber pressure, in head end of an SRM, decayed to 50 psia. 

6. Roll Maneuver — A roll maneuver of plus or minus 90 deg is required to place the vehicle in a 
heads-down or heads-up attitude, respectively, along the flight azimuth. The time to initiate the roll 
maneuver was fixed at 6.4 sec from launch, at a rate of 7.84 deg per second. These values were taken from 
Reference 8 and are consistent with the roll maneuver to heads-down currently used in STS trajectories. 
Note, in actuality the vehicle initiates the roll maneuver once a certain velocity is achieved. 

7 . Roll Maneuver to Heads-Down for an Initially Heads-Up Trajectory — The command to roll the 
vehicle to heads-down, for an initially heads-up trajectory, was fixed to initiate at a Mach number of 2.5, 
and a roll rate of 15 degrees per second. The reason for this roll maneuver is explained in the Program 
Modifications section. 

An example of all the namelist input required to run a trajectory using the MASTRE program is 
displayed in the Appendix. An explanation of all the variables listed is available in Reference 4. 


F. Procedure 

The procedure used in this study required numerous trajectory runs using the MASTRE program. 
First, SRM performance and mass property data were generated for sixteen pairs of motors, with propellant 
burnrates varying from 0.368 ips (current target) to 0.398 ips, at 0.002 ips increments. Note, these SRM 
data were generated from the information in the Groundrules Section 5, and each pair of SRMs had the 
same characteristics. Next, SRM data with 0.368 ips bumrate, no maximum Q constraint, and all the 
trajectory information from the Groundrules section were input into the MASTRE program. Then, the 
program was executed, and heads-up trajectory was optimized, based on maximizing injection weight, for 
this case. This process was repeated to attain optimal performance heads-up trajectories, with no max Q 
constraints, for SRM data with 0.370, 0.372, 0.374, 0.376, 0.378, 0.380, 0.382, 0.384, 0.386, 0.388, 
0.390, 0.392, 0.394, 0.396, and 0.398 ips. Next, the entire procedure is repeated implementing max Q 
constraints of 660, 680, 700, and 720 pounds force per square foot (psf). Therefore, a family of injection 
weight (payload capability) versus SRM bumrate curves are generated. Finally, for comparison, a heads- 
down trajectory using a 0.368 ips burnrate SRM and a 680 psf Q constraint was optimized for performance. 
This case incorporates the current flight conditions. 

Note, when maximum Q constraints were implemented, an SSME throttle down was required dur- 
ing first stage flight. For simplicity, the engines were allowed to operate at only 109 percent RPL and 65 
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percent RPL. Therefore, the time to initiate SSME throttle down to 65 percent RPL, for Q control, was 
based on satisfying the Q constraint while optimizing performance. Throttle up time (back to 109 percent 
RPL) occurred six seconds after max Q. 


III. RESULTS 


The results of this study are summarized in Figure 4, which is a plot of injection weight versus SRM 
bumrate. This figure shows the heads-up ascent, injection weight achievable for Q constraints of 660 psf, 
680 psf, 700 psf, 720 psf, and unconstrained Q; with SRM burnrates ranging from 0.368 ips to 0.398 ips by 
0.002 ips increments. The heads-down ascent (current flight mode), injection weight of 318,217 lbm, for 
an SRM bumrate of 0.368 ips, is shown for comparison. The unconstrained Q curve is essentially linear, 
with injection weight increasing with increasing SRM bumrate. The maximum injection weight value on 
this curve is 326,064 lbm, which occurs at an SRM bumrate of 0.398 ips and Q of 823 psf. For Q con- 
strained to 660 psf, the curve is nonlinear with injection weight increasing with SRM bumrate, to a maxi- 
mum value of 321,556 lbm at 0.390 ips. Thereafter, injection weight decreases with increasing SRM 
bumrate. The curve with a Q constraint of 680 psf is nonlinear, and injection weight increases with increas- 
ing SRM bumrate to a maximum value of 322,510 lbm, at 0.394 ips. Injection weight decreases with in- 
creasing SRM bumrate thereafter. The curves with Q constraints of 700 and 720 psf are also nonlinear. 
Both of these curves display injection weight increasing with increasing SRM bumrate, with respective 
maximum values of 323,416 lbm annd 324,103 lbm both occurring at 0.398 ips. 

Tables 1 , 2, 3, 4, and 5 contain the injection weight and SRM bumrate values used to generate 
Figure 4. Table 1 displays this information for unconstrained Q; whereas, Tables 2 through 5 display these 
data for Q constraints of 660 psf, 680 psf, 700 psf, and 720 psf, respectively. Tables 1 through 5 also list the 
following parameters: ignition weight, delta injection weight (heads-up value minus the baseline heads- 
down value), minimum product of dynamic pressure and angle of attack (Qa), minimum Qa time, max Q, 
max Q time, MECO time, and initiation of throttle down to 65 percent RPL time. Recall that the MECO and 
throttle down times are control parameters in the optimization scheme. 

The remainder of the results are mission elapsed time (MET) history plots of a variety of parameters 
and an altitude versus range plot. These figures compare the baseline heads-down ascent trajectory to the 
heads-up ascent trajectory with Q constrained to 680 psf and an SRM R b of 0.394 ips. This particular heads- 
up case was chosen for comparison to the baseline because 680 psf is a mission nominal design Q, and an 
SRM bumrate of 0.394 ips provided the optimum (maximum) injection weight at this Q. A discussion of 
each figure follows: 

Figure 5 — Compares the vacuum thrust-time traces for individual SRMs with burnrates of 0.368 
ips (used in the baseline heads-down trajectory) and 0.394 ips (used in the heads-up trajectory). The 0.394 
ips motor provides significantly more impulse prior to 80 sec than the 0.368 ips bumrate motor. In par- 
ticular, the total impulses at 60 sec are 170.034 x 10 6 lbf-sec and 186.500 x 10 6 lbf-sec for the 0.368 ips and 
0.394 ips motors, respectively. The difference between these values is 16.466 x 10 6 lbf-sec. However, at 
action time (motor burnout) the total impulses are 293.715 x 10 6 lbf-sec and 293.732 x 10 6 lbf-sec for the 
0.368 ips and 0.394 ips burnrates SRMs, respectively. These values differ by only 0.017 x 10 6 lbf-sec. 

Figure 6 — Shows SSME power level time histories during first stage, for the heads-up and heads- 
down trajectories. The second stage power level time history is constant at 109 percent RPL until initiation 
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Figure 4. Heads-up injection weight comparison. 



TABLE 1 . HEADS-UP PERFORMANCE/TRAJECTORY SUMMARY FOR UNCONSTRAINED Q 
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* BASELINE: HEADS-DOWN 
t RELATIVE TO BASELINE 

t NO THROTTLE DOWN REQUIRED FOR UNCONSTRAINED Q 






















TABLE 2. HEADS-UP PERFORMANCE/TRAJECTORY SUMMARY, Q CONSTRAINT OF 660 PSF 
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t RELATIVE TO BASELINE 

t THROTTLE LEVEL 100% RP' BFFOHE BUCKET, 109% AFTER BUCKET 




TABLE 3 . HEADS-UP PERFORMANCE/TRAJECTORY SUMMARY, Q CONSTRAINT OF 680 PSF 
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TABLE 4. HEADS-UP PERFORMANCE/TRAJECTORY SUMMARY, Q CONSTRAINT OF 700 PSF 
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TABLE 5. HEADS-UP PERFORMANCE/TRAJECTORY SUMMARY, Q CONSTRAINT OF 720 PSF 
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SRM VACUUM THRUST (LBF) 


0.368 IPS BURNRATE 

0.394 IPS BURNRATE 



Figure 5. SRM vacuum thrust comparison. 
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Figure 6. Heads-up versus heads-down SSME power level comparison. 
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of 3-G throttling. Recall, 100 percent RPL is roughly equivalent to a vacuum thrust of 469,760 lbf. The 
initiation of Q control throttle down to 65 percent occurs at 28.36 sec for the heads-down case, and at 20.33 
sec for the heads-up case. These times are optimized values, since the throttle down time is a control 
parameter in this study. Also, these throttle down times are consistent with the fact that the heads-up SRM 
provides more impulse early in flight, compared to the heads-down motor. 

Figure 7 — Is an altitude versus MET plot of the heads-down and heads-up trajectories. Of par- 
ticular interest is that the heads-up ascent flies a lofted profile relative to the heads-down ascent, until 
roughly 300 sec. After 300 sec the altitude profiles are essentially the same, because both trajectories have 
the same MECO altitude target (or terminal constraint). 

Figure 8 — Plots range versus MET for the heads-down and heads-up trajectories. Both range time 
histories have similar signatures. However, the heads-up ascent maintains a further down-range distance, 
relative to the heads-down ascent, from 100 sec to MECO. 

Figure 9 — Shows the altitude versus range flight profile of the heads-down and heads-up trajec- 
tories, thus combining Figures 6 and 7. Like Figure 6, examination of this plot shows the heads-up ascent 
flies a lofted profile relative to the heads-down ascent. 

Figure 10 — Displays the relative velocity (V R ) time histories of the heads-down and heads-up 
ascents. These relative velocity profiles are similar from SRB ignition to 50 sec, with the heads-up ascent 
velocity profile being consistently faster than the heads-down velocity profile. However, from 50 sec to 
approximately 120 sec, a significant difference exists between the velocity profiles. This difference is 
roughly 250 m/sec at 100 sec. The velocity profiles are again similar from roughly 120 sec to MECO; 
although, the heads-up trajectory maintains a higher velocity throughout this time period than the heads- 
down trajectory. However, the same target velocity is achieved at MECO in both cases. 

Figure 11 — Is a plot of relative flight path angle (y R ) versus MET of the heads-down and heads-up 
trajectories. This angle is measured in the pitch plane from the ET centerline to the local horizontal. 
Examination of the plot reveals that the heads-up ascent has larger values for y R until approximately 60 sec, 
relative to the heads-down ascent. From 60 sec to roughly 1 25 sec the heads-down ascent y R is greater than 
the heads-up ascent. Finally, from 125 sec to MECO the y R profiles are essentially the same because both 
trajectories must hit the same MECO target. Also, this plot is consistent with the heads-up ascent flying a 
more lofted trajectory than the heads-down ascent. 

Figure 12 — Shows the angle of attack (a) time histories of the heads-down and heads-up trajec- 
tories. This angle is measured in the pitch plane, from the ET centerline to the projection of the relative 
velocity vector. During first stage, after vertical rise and roll to flight azimuth maneuvers, these time 
histories should have the same signatures. The reason for this condition is due to the fact that in both trajec- 
tories, the same a-Mach number history is flown during the first stage pitchover (tiltover) sequence. 
Examination of the plot shows that the signatures are indeed the same from roughly 25 sec to 1 30 sec. From 
1 50 sec to 450 sec, the heads-up ascent flies a more negative a than the heads-down ascent. Then from 450 
sec to MECO, the heads-up ascent a is greater than that of the heads-down ascent. Note, the discontinuities 
in the a profile (where a jumps to zero degrees), between 100 and 130 sec, is due to an input problem. The 
a-Mach number history input table did not contain enough values required to get the vehicle to the staging 
time. Therefore, the MASTRE program assumed an a of zero once the input table ran out. However, all the 
a-Mach number history data available in Reference 8 was used in the input table. The problem would have 
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RANGE (NAUTICAL MILES) ALTITUDE (KILOMETERS) 



MISSION ELAPSEO TIME (SECONDS) 


Figure 7. Heads-up versus heads-down altitude comparison. 


DATA SOURCE: HDOWN368 RANGE HEADS-DOWN (SRM RB - 0.368 IPS) 



Figure 8. Heads-up versus heads-down range comparison. 
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RELATIVE VELOCITY (METERS/SECONDS) 


DATA SOURCE: HDOWN368 ALT HEADS-OOWN (SRM RB - 0.368 IPS) 
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Figure 9. Heads-up versus heads-down altitude/range comparison. 
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Figure 10. Heads-up versus heads-down relative velocity comparison. 
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RELATIVE FLIGHT PATH ANGLE (DEGREES) 


DATA SOURCE: HOOWN368 GAMR HEADS DOWN (SRM RB - 0.368 IPS) 

DATA SOURCE: HEADUP394 GAMR HEADS-UP (SRM RB ■ 0.394 IPS) 
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Figure 1 1. Heads-up versus heads-down relative gamma comparison. 
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been alleviated by adding additional Mach number values to the table, which holds the last a value in 
Reference 8. Obviously, these discontinuities will be observed in the Qa time history, and will also 
manifest in the \p time history to a lesser degree. Since this input problem is of short duration and does not 
effect the min-H attitude control phase of the program, no significant difference would be observed in the 
optimized injection weight. Therefore, no updated trajectories were generated. 

Figure 1 3 — Displays the angle of sideslip ((3) versus MET for the heads-down and heads-up trajec- 
tories. Note, 3 is measured from the ET centerline to the projection of the relative velocity vector, in the 
yaw plane. The two 3 profiles, like the a profiles, have the same signatures during the pitchover phase from 
roughly 25 sec to 1 30 sec. This condition results because the same (3-Mach number history is flown during 
first stage pitchover sequence. The second stage profiles, from approximately 130 to MECO, are also 
similar. However, the heads-up ascent maintains greater values of (3 relative to the heads-down ascent, 
during this time period. Note, discontinuities (where (3 jumps to a value of zero) exist in the (3 time his- 
tories. This situation, like the a time history discontinuities, is due to the lack of sufficient data being 
supplied in the input (3-Mach number history table. These discontinuities will obviously be present in the 
Q(3 time history and will also manifest itself, to a smaller degree, in xy- Again, since this problem is of 
short duration and prior to the min-H attitude control phase, no significant difference would be observed in 
the optimized injection weight. Therefore, no updated trajectories were generated. 

Figure 14 — Shows the Mach number time histories of the heads-down and heads-up trajectories. 
The Mach number profiles have similar signatures as their respective V R profiles, as expected. The Mach 
number histories are similar from SRB ignition to 50 sec, with the heads-up ascent having consistently 
higher Mach number values relative to the heads-down ascent. However, significant differences exist 
between the two profiles from 50 sec to roughly 1 20 sec. This difference is about 0.50 at 1 00 sec. The Mach 
number profiles are again similar from approximately 1 20 sec to MECO; with the heads-up trajectory main- 
taining higher Mach number values relative to the heads-down trajectory during this time frame. 

Figure 15 — Is a plot of pitch attitude angle (xp) versus MET for the heads-down and heads-up 
trajectories. The heads-up trajectory has substantially smaller values of Xp> relative to the heads-down 
trajectory, from SRB ignition to roughly 50 sec. This situation is consistent witth lofted trajectory of the 
heads-up ascent relative to the heads-down ascent. From 50 sec to roughly 1 25 sec, the heads-up ascent has 
larger xp values than the heads-down ascent. After 125 sec, both Xp time histories are essentially the same. 
Recall, Xp is a control parameter in the MASTRE program and therefore, these profiles are optimized to 
maximize vehicle injection weight. 

Figure 1 6 — Shows the yaw attitude angle (xy) time histories for the heads-down and heads-up 
ascents. The xy profiles are significantly different from SRB ignition until about 20 sec, due to one profile 
rolling to heads-down while the other rolls to heads-up. Note that roll and yaw are coupled. The profiles are 
basically similar from 20 sec to approximately 125 sec, with the heads-up ascent having consistently larger 
values of xy compared to the heads-down ascent. However, at roughly 80 sec, the heads-up xy time history 
has a local peak which does not appear in heads-down ascent. This peak occurs because during this time 
frame the vehicle achieves a Mach number of 2.5 and thus triggers the heads-up ascent to roll to heads- 
down. The reason for the roll to heads-down, for an initially heads-up trajectory, is given in the Program 
Modifications section. Finally, from roughly 125 sec to MECO, the Xy values for both cases are essentially 
the same. Also, recall xy is optimized by the MASTRE program to maximize vehicle injection weight, 
since it is a control parameter. 
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ANGLE OF SIDESLIP (DEGREES) 



MISSION ELAPSED TIME (SECONDS) 


Figure 13. Heads-up versus heads-down beta comparison. 



MISSION ELAPSED TIME (SECONDS) 

Figure 14. Heads-up versus heads-down Mach number comparison. 
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Figure 17 — Displays the roll attitude angle (\r) time histories for the heads-down and heads-up 
trajectories. The roll maneuvers, after launch pad clearance, can readily be observed in the plot from 6.4 
sec to roughly 18 sec for both ascent profiles. The heads-down ascent rolls from \r of 90 deg to 180 deg; 
whereas, the heads-up ascent rolls from 90 deg to 0 deg. Also, for the heads-up ascent, the roll maneuver to 
heads-down can be seen in the plot from roughly 80 sec to 92 sec, where Xr ramps from 0 deg to 1 80 deg. 
Recall, this second roll maneuver increases performance for reasons stated in the Program Modifications 
section, and is triggered to initiate when the vehicle attains a Mach number of 2.5 . Note, the discontinuities 
in the plot, where the value of xr jumps from 180 deg to -180 deg (or vice-versa), are misleading. A xr 
value of 180 deg or -180 deg, after the roll maneuver to flight azimuth, are actually the same attitude 
position. 

Figure 18 — Is a plot of dynamic pressure (Q) versus MET for the heads-down and heads-up 
ascents. Only the first stage time history of Q is shown because Q is essentially zero during second stage. 
Both ascent cases have the typical bell shaped Q profile, and they were both constrained to the same 
maximum value of 680 psf. However, max Q occurs at 45.75 sec for the heads-up trajectory and at 52.04 
sec for the heads-down trajectory. Therefore, a difference of 6.29 sec is observed in the max Q time values 
between the two cases. 

Figure 1 9 — Shows the product of dynamic pressure with angle of attack (Qa) time histories for the 
heads-down and heads-up trajectories. Only the first stage Qa time histories are shown because Qa is 
essentially zero during second stage. The signatures of these Qa profiles are the same because both ascent 
cases are required to fly the same a profile during the tiltover sequence. Note, a tiltover control initiates at 
18.4 sec and continues through staging; whereas, initial tiltover from 6.4 sec to 18.3 sec, is optimized by 
the MASTRE program. Also, the minimum Qa values of -3049 psf-deg and -3052 psf-deg for the heads- 
down and heads-up ascents, respectively, are approximately the same. However, the heads-up ascent 
minimum Qa occurs 6.5 sec earlier than that of the heads-down ascent. These minimum Qa values are 
consistent with the PRM-4 design Qa of -3000 psf-deg, as stated in Reference 7. The discontinuities in the 
Qa time histories (where Qa jumps to a value of zero psf-deg), between 100 and 130 sec, is due to the input 
a-Mach number history table not being of sufficient length. 

Figure 20 — Displays the product of dynamic pressure with angle of sideslip (Q(3) versus MET for 
the heads-down and heads-up trajectories. Only the first stage values of Q(3 are shown because Q(3 is essen- 
tially zero during second stage. Both trajectories have similar Q(3 profiles during the tiltover sequence, 
from 1 8 . 3 sec to staging , because the same (3-Mach number history is input . However, during the initial rol 1 
maneuvers, from 6.4 sec to approximately 18 sec, large differences exist in the Q(3 profiles. The heads- 
down ascent achieves a minimum Q(3 value of roughly -1250 psf-deg; whereas, the heads-up only drops to 
a minimum of approximately -250 psf-deg. This difference occurs because one trajectory case rolls to 
heads-down and the other rolls to heads-up, and yaw and roll are coupled. The discontinuities in the Q(3 
profiles (where Q(3 jumps to 0 psf-deg), between 100 and 130 sec, is due to insufficient length of the input 
(3-Mach number history table. 

Figure 21 — Is a plot of axial (longitudinal) acceleration versus MET for the heads-down and 
heads-up ascents. The acceleration profiles during second stage, from roughly 125 sec to MECO, are 
almost identical. Note, the maximum vehicle acceleration was constrained to 3 G’s in both cases. How- 
ever, the first stage time histories are significantly different, with the heads-up ascent having higher values 
of acceleration than those of the heads-down ascent. This condition is a direct result of higher thrust pro- 
vided by the 0.394 ips bumrate SRM relative to the 0.368 ips bumrate SRM (Fig. 5). 
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ROLL ATTITUDE ANGLE (DEGREES) 








OATA SOURCE: HDOWN368 QALPH HEADS-DOWN (SRM RB . 0.368 IPS) 



Figure 19. Heads-up versus heads-down Q*alpha comparison. 



MISSION ELAPSED TIME (SECONDS) 


Figure 20. Heads-up versus heads-down Q*beta comparison. 
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DATA SOURCE: HDOWN368 
DATA SOURCE: HEADUP394 


— AXACC 
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HEADS-DOWN (SRM RB • 0.368 IPS) 
HEADS-UP (SRM RB . 0.394 IPS) 

i i i ■ f 



Figure 21. Heads-up versus heads-down acceleration comparison. 


IV. DISCUSSION 


The evolutionary path of the MASTRE program began with the ROBOT program developed by Dr. 
R. G. Gottlieb [2]. ROBOT was used to optimize the upper stages of the Saturn rocket series, and since 
these rockets had symmetric features, a near optimum attitude profile could be attained by the use of a 
gravity (zero angle of attack) profile. Nonatmospheric optimization algorithms could be used on Saturn 
upper stages, because its trajectory was significantly out of the atmosphere at both booster cutoff and orbit 
insertion. However, the airplane-like design and low injection altitude of the Space Shuttle requires the 
calculation of aerodynamic, atmospheric, asymmetric center of gravity, and thrust effects during the 
complete ascent flight trajectory. In other words, the STS complexity and trajectory requirements necessi- 
tated implementation of atmospheric optimization, along with the additional constraint of thrust vector 
control moment balance. The Lifting ROBOT program [3], also developed by Dr. Gottlieb, and the 
RAGMOP program [1], developed by a team assisted by Dr. Gottlieb and Dr. J. T. Lyons, were the first 
optimization attempts within the atmosphere and STS related constraint enforcement. Features from these 
programs are incorporated into the current MASTRE program, developed by a team led by Dr. Lyons, 
along with the additional adjoint equations required by the complex equations of motion of the STS. 

The min-H optimization strategy employed by the MASTRE program has proved to be a computa- 
tionally fast and stable method of determining optimum trajectories. Most simulation cases converged to 
solutions in seven to twelve iterations; however, convergence problems existed for heads-up trajectories 
initially. The MASTRE program was written to allow the Space Shuttle to fly in either a heads-down or 
heads-up mode; but since the heads-down mode is the current flight design, only this option had been vastly 
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exercised. Essentially, this study is the first time the heads-up option in the MASTRE program required 
use. Therefore, some minor problems were encountered when first attempting to fly the Shuttle heads-up. 
Some of these problems included: ( 1 ) verifying that the heads-up option flag was indeed causing a heads-up 
attitude, (2) determining the correct sign on roll rate so as to cause a 90 deg roll to heads-up instead of a 270 
deg roll, (3) determining the cause of errors in the integration subroutine, and (4) selecting initial 
estimates of the control parameters which would allow the program to converge to a solution. Once these 
problems were solved, a heads-up attitude could be flown from the start of the pitchover phase through 
MECO. However, recall that a 180 deg roll to heads-down prior to SRB jettison, initiated at Mach 2.5, 
improved the performance of an initially heads-up trajectory. So, a new code was added to the MASTRE 
program to allow for this second roll maneuver if desired. This modification resulted in the same latter two 
minor problems stated earlier for an entire heads-up ascent. Note that unless otherwise stated, a heads-up 
ascent in this study refers to a heads-up attitude from initiation of tiltover until the vehicle achieves a Mach 
number of 2.5, where a roll maneuver to heads-down is initiated and is maintained thereafter. The problems 
associated with the additional roll maneuver were eventually solved once proper initial control parameter 
values were selected. Therefore, the heads-up trajectories were able to converge to an optimum solution. 

The results of this study indicate substantial improvements in the STS lift capability can be achieved 
if heads-up flight is combined with increased SRM ballistic performance. One of the major drivers regulat- 
ing the amount of increase of STS performance is the maximum Q limit. If Q is unconstrained, then STS lift 
capability is limited only by propulsion system performance; neglecting vehicle structural, load, and 
thermal limits. For example, a heads-up trajectory, with 0.398 ips bumrate SRMs and max-Q of 823 psf, 
increases injection weight by 7847 lbm compared to the baseline heads-down trajectory, with 0.368 ips 
bumrate SRMs and 680 psf max-Q. However, the dispersed orbiter thermal protection system (TPS) Q 
limit is 8 19 psf. Therefore, mission design Q’s are nominally around 680 psf in order to provide protection 
against three sigma system and wind dispersions. So, with Q constrained to this maximum design value, an 
optimum SRM bumrate of 0.394 ips results for heads-up trajectories. Injection weight is increased by 4293 
lbm under these conditions compared to the baseline heads-down ascent. Note that the total injection 
weight values listed in Tables 1 through 5 do not accurately reflect potential STS performance capability as 
of this writing. However, the delta injection weight values would be very similar. One reason for the inac- 
curacy is because an SSME power level of 109 percentRPL was used in all trajectory runs, and this level is 
currently not to be used except during an abort mode to avoid orbiter ditching. A great deal of testing and 
analysis needs to be conducted before a 109 percent SSME power level setting is certified for nominal 
ascent. The other reason for inaccurate injection weights results from the use of FWC SRM’s on all trajec- 
tories. This type of SRM was under development and testing when this study was initiated. However, due 
to the STS 5 1 -L (Challenger) accident and the ensuing SRM redesign, development of this motor has been 
mothballed. Also,, studies are now underway to examine the use of an Advanced SRM (ASRM) on the 
STS. Additional comments on the ASRM are mentioned later in this discussion. 

A heads-up trajectory with 0.368 ips bumrate SRMs, does not offer a performance advantage 
compared to the baseline heads-down trajectory. Therefore, payload lift capability of the STS is unaffected 
for heads-up ascent with no increase in SRM ballistic performance. However, SSME throttling is no longer 
necessary under these circumstances to satisfy a max-Q constraint of 680 psf. In fact, max-Q is only 662 psf 
for the heads-up ascent with 0.368 ips bumrate SRMs and unconstrained Q. So, during first stage the poten- 
tial problems associated with SSME hydraulic lock-up at low power levels are eliminated. Note, this trajec- 
tory’s reduced Q value of 662 psf does not produce significant load advantages, because the Qa profile is 
essentially unchanged from the baseline heads-down trajectory. 
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The performance gain achieved with heads-up flight is due to the relationship between optimum 
SRM bumrate and Q and Qa constraints being modified to accommodate higher SRM bumrates. The 
change in the relationship occurs because the existing Qa constraint forces heads-up ascents to fly at a 
higher altitude or become more lofted, relative to a heads-down ascent (Figs. 7 and 9). The higher altitude 
profile of the heads-up trajectory decreases dynamic pressure, which decreases SSME throttling 
requirements and allows flight with higher performing SRMs. These SRMs provide the increased impulse 
necessary to cause a heads-up ascent to be faster, as well as higher, at SRB separation than a heads-down 
ascent. Therefore, the overall vehicle performance increases due to increased first stage performance. 

Although heads-up flight with increased SRM propellant bumrate enhances STS performance, it 
also poses several potential problems. One of the most significant problems concerns lift-off load 
increases. Maximum SRM pressure rise rate and maximum SRM thrust rise rate increase with faster 
propellant bumrates where the larger rise rates result in increased lift-off loads. Heads-up feasibility studies 
conducted by Rockwell International [6] indicate significant exceedances in lift-off loads for 0.392 ips 
bumrate motors. This study also revealed slight exceedances in the post high-Q flight regime loads for 
0.392 ips bumrate SRMs. Other issues associated with heads-up flight involve initialization load (I-load) 
updates and onboard computer software additions. According to Reference 6, a detailed flex/rigid body 
stability analysis is required to determine new I-loads, such that control stability margin is maintained. 
Also, to accomplish a roll maneuver to heads-down at Mach 2.5, a roll table needs to be added to the 
onboard general purpose computer (GPC) flight software, along with elevon schedule changes. A final area 
of concern deals with ascent venting. Compartments in the STS are designed to vent with increasing alti- 
tude in order to equalize pressure. The lofted trajectory associated with heads-up flight may require 
compartment venting to occur faster than the capabilities of the current system. Further analysis is required 
to determine if ascent venting is indeed a problem. Overall, some of the heads-up issues mentioned are 
minor and can be resolved; however, the lift-off load exceedance problem indicates SRM propellant 
bumrate cannot drastically change. 

SRM ballistic performance can be increased by several methods besides increasing propellant 
bumrate. These methods include increasing propellant load, using a different propellant, and/or changing 
internal grain geometry. Assessment of these options is currently being evaluated under ASRM contracts 
managed by MSFC. Note, of primary importance in the ASRM studies is to determine methods to increase 
safety, reliability, and reproducibility of the SRM. Finally, recent studies by Rockwell International [13] 
show that if the SRM thrust-time trace shape is recontoured and bumrate is maintained at 0.368 ips, a 
heads-down ascent can improve the STS performance more than the heads-up ascent. Therefore, the heads- 
up flight mode does not always offer a performance advantage compared to the heads-down mode. 

The Space Shuttle is a complex vehicle and any change in one area usually affects many other areas. 
This study examined potential STS performance enhancements due to heads-up flight, and showed 
performance improvements are possible. However, the method used to increase STS performance was to 
increase SRM bumrate, which adversely affected lift-off loads. Further studies are being conducted to 
determine a boost phase propulsion system design, which improves STS performance without causing 
significant side-effects to the overall system. Whether the vehicle flys in a heads-up mode or continues to 
fly in a heads-down mode may well depend on the booster propulsion system it employs. 
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V. CONCLUSIONS 


From the results of this study, the following conclusions can be drawn: 

1. Heads-up and heads-down STS ascents can be flown with the same a-Mach number and 
(3-Mach number histories. 

2. The optimum SRM propellant bumrate for STS heads-up ascent, with a mission nominal design 
Q constraint value of 680 psf, is 0.394 ips. 

3. Heads-up ascent, with optimum burnrate SRMs and a 680 psf Q constraint, enhances perform- 
ance by 4293 lbm compared to the baseline heads-down ascent, with 0.368 ips bumrate SRMs and a 680 psf 
Q constraint. 

4. Heads-up ascent employing current target, 0.368 ips, bumrate SRMs does not offer a perform- 
ance advantage relative to the baseline heads-down ascent. However, SSME throttling for Q control is not 
necessary in this case, which is advantageous because during first stage the risk of engine hydraulic lock-up 
at low power levels is eliminated. 
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